Magnetic anisotropy in ferromagnetic Josephson junctions 
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Magnetotransport measurements were done on Nb/A^Os/Cu/Ni/Nb superconductor- insulator- 
ferromagnet-superconductor Josephson tunnel junctions. Depending on ferromagnetic Ni interlayer 
thickness and geometry the standard (lei) magnetic field dependence of critical current deviates 
from the text-book model for Josephson junctions. The results are qualitatively explained by a 
short Josephson junction model based on anisotropy and 2d remanent magnetization. 
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Superconductivity (S) and ferromagnetism (F) in thin 
layered films have now been studied during some decades 
[ ]. In SF bilayers the superconductivity may be non- 
uniform [2], i.e. the Cooper pair wave function extends to 
the ferromagnet with an oscillatory behavior. In Joseph- 
son junctions (JJs) based on s-wave superconductors the 
phase coupling between the superconducting electrodes 
can be shifted by 7r when using a ferromagnetic barrier 
with an appropriate chosen thickness dp, i.e. SFS or 
SIFS-type junctions (I: insulating tunnel barrier). Only 
in recent years the experimental realization of tt JJs was 
successful. In particular, the tt coupling was demon- 
strated by varying the temperature [3, 4, 5], the thickness 
of the F-layer [5, 6, 7, 8] or measuring the current-phase 
relation of JJs incorporated into a superconducting loop 
[9, 10, 11]. The coupling can also change within a single 
JJ by a step-like F-layer, i.e. one half is a JJ and the 
other half is a tt J J [12, 13]. 

For useful classical or quantum circuits based on 
SFS/SIFS JJs a large critical current density j c (small 
Josephson penetration depth Aj) and a high I C R prod- 
uct are needed [14, 15]. Up to now the limiting fac- 
tor is the low jc due to strong Cooper pair breaking in- 
side F-layer. Alloys of magnetic and non-magnetic atoms 
such as NiCu face problems of clustering [16] and strong 
magnetic scattering [5, 8]. Promising experiments using 
strong ferromagnet transition metals [6, 17, 18, 19, 20] 
were published. 

Shape anisotropy of magnetic interlayer may provoke a 
not flux-closed domain structure and consequently a shift 
of critical current diffraction pattern I C {H) [19]. In ex- 
periments [3, 5, 7, 19, 20] the SFS/SIFS JJs had nearly 
mirror-symmetrical I C (H), i.e. the effective shift along 
H-axis is small, usually less than one flux quantum $0. 
This is explained by a multi-domain state of F-layer with 
a very small net magnetization. However, up to now the 
2d nature of thin-film magnetism was disregarded. 
In this Letter the I C (H) dependence for remanent 2c? 
magnetization of F-layer is systematically studied. First, 
the maximal flux from F-layer is estimated. Second, the 
Ic(H) pattern considering 2d in-plane magnetization is 
calculated for different aspect ratios. Third, the I C (H) 
pattern is measured along both magnetic axes for various 
junction geometries and dp. 
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FIG. 1: (Color online) Calculated surface plot of I C (H X , H y ) 
for different geometries and remanent magnetizations a. The 
magnetization vector M points from top right to bottom left 
corner (arrow), 7° is shifted in opposite direction. 



The maximal shift of I c (H y ) is estimated for a strong 
magnet, i.e. Ni, being magnetized fully in-plane and 
along j/-axis. The atomic magnetic momentum is 0.6 /ib 
[21], the specific density p is 8.9 g/cm 3 (bulk) and mag- 
netization /itjM = 0.64 T. A cross-section of length 
L x = 100 /iin and F-layer thickness dp = 3 nm encloses 
a magnetic flux $m = dpL x pu^M . The total magnetic 
flux $ through the JJ is the applied field flux = 
{2\i J + dp)L x H y (London penetration depth A^ = 90nm) 
plus $ A f, i.e. $ = $ H ± $ M = 8.85$o#y/mT ± Q2$ . 
The I c (H y ) pattern is shifted by 92 periods from the 
center. This simple calculation neglects dead magnetic 
layer [22], as found in SFS/SIFS JJs [5, 8, 20], and de- 
magnetizing by domains. Real Ni films tend to form com- 
plex magnetization profiles (in/out-of-plane) and domain 
structures as function of dp [23]. Integral magnetization 
measurements in SF layers show a complex behavior as a 
function of temperature, applied field, and sample history 
[24, 25], e.g. SF structures spontaneously alter their stray 
field by changing magnetic domain distribution [26]. The 
local magnetization depends on stray fields from neigh- 
bor domains, flux focusing from S-electrodes [19, 27] or 
on bias induced spin accumulation at F/S interface [28]. 
It is generally agreed that the average magnetization in 
SFS/SIFS JJs is much smaller than the maximum mag- 
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FIG. 2: (Color online) Calculated I c as function of Id mag- 
netic field, I c (H x ,0) (top) and I c (0,H y ) (bottom) for the 
three geometries and various a. A substantial deviation from 
ideal (a = T/m) pattern appears for already small magne- 
tization vector M = a(L x , L y ). 
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FIG. 3: (Color online) Measured I c (H y ) pattern of SIFS JJs 
for thin (top) and thick (bottom) Ni layer and different ge- 
ometries. Ic(Hx) pattern of thick Ni layer is plotted in gray. 
Onset of magnetic anisotropy effects in Ni layer is between 
2.2 and 3.9 nm. 



netization estimated above. But as shown below even a 
remanent 2d magnetization of 1% of the maximal value 
may notable change the I C (H) pattern. 

A qualitative model for I C (H) in presence of a uniform, 
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M = a(L x ,L y ) is assumed to be orientated in-plane 
(Meissner-screening of S-electrodes) along the diagonal 
(i.e. longest) axis of sample (rough approximation for 
magnetic shape anisotropy). The magnetization \M\ is 
several orders of magnitude smaller than the upper limit 
given by a fully saturated magnetic layer. Note that the 
easy axis of ferromagnetic film can be determined by the 
magnetic field during deposition, too. The model for <t>jv/ 
is just exemplary for the effect of 2c? in-plane magnetiza- 
tion in I C (H). 

In Fig. 1 the surface plot of I c (H x ,H y ) is depicted for 
various a and geometries. The position of I® is shifted 
from the center (H x = H y = mT) in opposite direction 
of M. Fig. 2 depicts I c (H x ,0) and I c (0,H y ) pattern. 
These graphs resemble standard Id I C (H) (H = H x ,H y ) 
measurements. For some \M\ / 0T a single-peaked 
I C (H X ), I c (H y ) pattern is calculated, which -on first 
glance- resembles a shifted \sin(H)/ H\ Fraunhofer pat- 
tern. However, the maximum I c is smaller than the 
real I® and the height of side-maxima do not obey the 
expected value. For example the I c (H x ,0) pattern of 
30 x 30 fim 2 sample and a = 0.071 T/m is shifted by 
less than <t>o and its maximum I c is already reduced to 
~ 0.8/?. This simple I C (H X , H y ) model may qualitatively 
explain some experimental observations (Fig. 3) on SIFS- 
type JJs with Ni interlayer. 

For experiment JJs with similar areas, i.e. 30 x 30, 50 x 
100 and 25 x 200 /im 2 , were fabricated. The deposition 
and structuring of JJs is described in Ref. 29. The 
SIFS multilayer was magnetron sputter deposited with 
Ni thickness dp ranging from 1—6 nm. The tunnel bar- 
rier was formed for 30 min at a partial oxygen pressure 
of 0.1 mbar. After oxidation a 2 nm Cu film was in- 
serted. All JJs were deposited in a single run by shift- 
ing the substrate and target to obtain a wedge-shaped 
Ni-layer. Normal state and subgap resistance indicate a 
small junction to junction variation. The IV and I C (H X ), 
I c (H y ) characteristics were measured at 4.2 K for two 
sets of samples {dp = 2.2 , 3.9 nm). Cooldown was done 
in zero field and thermal cycling up to w 15 and 300 K 
to check reproducibility. Transport measurements were 
made in a liquid He dip probe using low-noise home made 
electronics and room-temperature voltage amplifier. The 
magnetic fields (H x ,0), (0,H y ) were applied in-plane of 
the sample and parallel the junctions axis (Fig. 1). The 
voltage criteria V c for I C (H) determination was 0.2— 1/iV. 
A lower subgap resistance for dp = 3.9 nm sample leads 
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to larger offset currents. Positive and negative current 
branch of 1VC had similar magnetic field dependence 
+I c (H y ) ~ | — I c (H y )\. Magnetic field was swept between 
±1.5 mT for all samples. All junctions had their lateral 
sizes comparable or smaller than the Josephson penetra- 
tion length Aj, except the longest sample (dp = 2.2 nm, 
25 x 200 /im 2 ), whose L y is not strictly inside the short 
J J limit. 

The <1f = 2.2 nm samples showed very regular I c (H y ) 
pattern. All maximum 7 c 's were nearly centered and 
the spread of j c was ~ 1.7%, as determined from the 
maximum J c 's. The I C (H X ) pattern is symmetric, too 
(not shown). The oscillation period of I c {H y ) were de- 
termined by magnetic cross-section ~ 1/L X , and nearly 
independent of aspect ratio. No indication for a distorted 
supcrcurrcnt transport due to alloying at the Nb/Ni 
interface [30] can be found. Effects due to magnetic 
anisotropy were not detectable, either because the sam- 
ples were still inside the dead magnetic layer, or the 
anisotropy was absent or totally out-of-plane. 
The dp = 3.9nm samples had completely different I c (H y ) 
pattern showing in-plane magnetic anisotropy with some 
characteristic features. All maximum 7 c 's were shifted 
from the center, and the amplitude of shift increased 
with L y , i.e. « 0.24 mT for 30 x 30 /im 2 , w 0.5 mT for 
50 x 100/mr 2 , and w 0.8mT for 25x200/im 2 samples. The 
direction of shift varies between samples even if they were 
cooled and measured at the same time (random polarity 
of magnetic configuration). The position of main peak of 
I c (H y ) was reproducible after thermal cycling to 300 K. 
The width of main maxima (measured at offset line) was 
not strictly ~ 1/L X , and varies from sample to sample. 
The pattern were asymmetric, i.e. the height of same- 
order side maxima differed, probably due to non-uniform 
flux guidance in F-layer and re-orientation of domains. 
By rotating the magnetic field by 90°, i.e. measuring in 
Ic[H x ) mode, low 7 c 's, being nearly independent of H x , 
were detected. Even the squared shaped 30 x 30 /im 2 
J J had an almost flat I C {H X ) pattern. This indicates 
some magnetic crystallographic anisotropy along y-axis, 
probably caused by magnetron sputter deposition. Small 
deviations of 7 C (0, 0) for I C (H X , 0) and 7 C (0, H y ) measure- 



ments can be related to variations of magnetic configu- 
ration by the unshielded sample handling at 300 K. A 
considerable spread of maximum I c (Fig. 3) can be al- 
ready seen for JJs with same geometry, which is even 
increased by considering the maximum j c for different 
geometries. Simulations (Fig. 2) show that already a 
moderate magnetization M (a < 0.1 T/m) yields very 
different maximum 7 c 's of I C (H X ) and I c (H y ). A sample 
to sample variation of direction and amplitude of M ex- 
plains the data spread of Fig. 3. Obviously, this leads to 
very large variations in the I c {dp) dependence. 
In literature, the I C (H) pattern of SFS/SIFS JJs with 
comparable strong magnets were either shown for sam- 
ples with thin dp [19], or had deviations from the ideal 
sin(77)/77| form. For example, the maximum I c in Ref. 
20, Fig. 5 (inset) is too small compared to the first side- 
maxima. These samples were small in area (^ 1 /^m 2 ), 
and the F-layer could have been in single domain state. 
As 7 C varied smoothly with H, either M rotated softly, 
or a multi-domain structure with averaged 2d magneti- 
zation existed. For both cases the 2d nature of remanent 
magnetization may have suppressed the maximum I c - 
determined from I C (H) pattern- below the real 7°. 

In summary, the I C (H) pattern along both field axis 
of SIFS JJs with Ni interlayer were measured. Assum- 
ing magnetic anisotropy the characteristic features, i.e. 
shift or absence of central peak, can be qualitatively re- 
produced by simulations. As conclusion, the Id I C {H) 
pattern in presence of magnetic anisotropy can not yield 
the real 7°. Future experiments on SFS/SIFS JJs should 
be done by 2-axis I c (H x ,H y ) scan of JJs with well- 
controlled magnetic interlayer domain configuration. Su- 
perconducting spintronic devices such as FSF spin valves 
and 7r, 0-7T JJs are promising candidates for future clas- 
sical and quantum computers. As shown by this Letter 
control of magnetic anisotropy and magnetic domain con- 
figuration is essential for phase-switchable FSF/SFIFS or 
SFS/SIFS devices. 
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